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STFOPSIS 

The present work 'deals with, static flow-shop schedul- 
ing problem. Many exact and heuristic procedures have been 
reported in the literature for this purpose considering 
minimization of make span time as the optimization criterion. 
The exact procedures yield optimal solution but tend to be 
highly computationally inefficient for large sized problems 
involving more than two machines. The reported heuristic 
procedures, though computationally fairly efficient, do 
not guarantee the optimality of the solution. 

Given a sub optimal solution obtained by using one of 
the known heuristics, one can possibility , exploit the 
characteristics of the solution for developing heuristic 
rules which when used on this sub-optimal sequence would 
yield improved solution. In a broad sense, this is the 
basic idea used in the present work. 

The given suboptimal solution would result in idle 
time for some jobs as well as the machines. The information 
about the idle time has been to an advantage for the deve- 
lopment of proposed heuristic which is used for improving 
upon the suboptimal .sequence. The concepts of Earliest 
and Latest Finish Gantt matrices, slack time and critical 
path have been introduced for the development of the 
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heuristic. The heuristic is iterative in nature and in 
each iteration it identifies the jobs which become the 
candidates for the cyclic pairwise exchange. The iterative 
procedure is continued till it fails to generate any fur- 
ther improved sequence. 

Computer programmes for the Petrov's algorithm to 
obtain the initial sequence and for the proposed heuristic 
procedure have been developed. The performance of the 
proposed heuristic has been tested on a number of randomly 
generated problems of various sizes. It was observed that 
the heuristic improved the initial sequence for about 
70 percent of the test problems, further on an average 
there was a reduction of about 8 percent in the make span 
time for the problems in which the proposed heuristic was 
successful in generating the improved sequence. 



CHAPTER I 


INTRODUCTION 


Scheduling is the allocation’ of resources over 
time to perform a collection of tasks. This rather 
general definition when conceived in the context of 
manufacturing, assumes the resources to be machines 
and the tasks to be jobs or products that have to be 
processed. 

The term scheduling conveys two different mean- 
ings that are important to be mentioned here. First, 
scheduling is a decision-making function; it is the 
process of determining a schedule. Second, scheduling 
is a body of theory; it is a collection of principles, 
models, techniques and logical conclusions that provide 
insight into the scheduling function. 

The other term that is used quite frequently 
in this context is ' sequencing' . Sequencing is concerned 
with the arrangements and permutations in which a set 
of jobs under consideration are performed on all machines. 
Scheduling is concerned with the specification of start- 
ing or completion times of jobs on all machines. Sequenc- 
ing decisions focus on the arrangements of events, whereas 
scheduling decisions focus on the time of events. 



The present work deals with, the flow-shop schedul- 
ing problem. In a flow-shop, jobs are multistage in 
nature and the machines have serial configuration. A 
job is considered to be a collection of operations in 
which a special precedence structure applies. In parti- 
cular each operation after the first has exactly one 
direct predecessor and each operation before the last 
has exactly one direct successor. This is depicted in 
fig. 1.1 for a job j. The numbers 1,2,..., m, refer to 
the machines and each circle denotes an operation. Each 
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operation requires a different machine. A flow-shop is 
characterized by unidirectional flow of work, or in 
other words, it is possible to number the machines so 
that if the i-th operation of any job j precedes its 
k-th operation then the machine required by the i-th 
operation has a lower number than the machine required 
by the k-th operation. The machines in a flow-shop 
are numbered 1 , 2 ,..., m and the operations of job j 
are correspondingly numbered (l,j), ( 2, j ),..., (m,j). 
Each job can be treated as if it had exactly & opera- 
tions, for in cases where fewer operations exist, the 
corresponding processing times can be taken to be zero. 
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Given a flow-shop environment with m machines 
and n jobs, the objective is to deteimine an optimal 
sequence considering some measure of effectiveness. The 
effectiveness of a sequence can be measured in terms 
of makespan time (total elapsed time), average comple- 
tion time, due date performance, machine utilization, 
inventory of jobs in process, etc. Considerable research 
effort has been directed by many researchers towards the 
solution of above stated combinatorial optimization 
problem under a set of standard assumptions (l). 

Most of the researchers have considered the 
minimization of makespan time as the criteria, for gene- 
ration of optimal sequence. For a 2-machine, n job 
situation Johnson has presented an algorithm for 
obtaining an optimal solution. As the number of machines 
increases beyond two, the use of operation research tech- 
nique which can yield optimal solutions becames restric- 
tive due to large computational requirements. Therefore, 
number of heuristic procedures have been suggested which 
though computationally efficient can only yield near 
optimal solutions to the problem. The near optimal 
sequence obtained by any one of the available heuristic 
will, in general, have some idle times associated with 
its operations. The immediate question which comes to 
ones mind is, can these idle-times be exploited for the 
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generation of a tetter sequence? The present work 
essentially concerns itself with the generation of 
tetter sequence through the use of idle-time informa- 
tion. The criteria of minimizing make-span time is 
used for the evaluation of the sequences. Further, 
the standard set of assumptions listed by Bakshi and 
Arora (1) are considered for the present work. 


The general frame work of the proposed proce- 
dure is depicted in Fig. 1.2. 



Fig. 1.2 


This is a kind of feed-back system in which the initial 
sequence can be any sequence. The algorithm attempts 
to improve upon the given sequence S. Each application 
of the algorithm yields a new sequence. If the new 
sequence S* happens to be better than the previous one, 
the algorithm is applied again on the new sequence. The 
procedure continues till no further improvement is 
possible. 

It should be noted here that better the initial 
sequence, faster is the convergence towards final 
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improved sequence for the problem under consideration. 
Therefore, the initial sequence is generated by apply- 
ing any one of the available heuristics instead of 
choosing it randomly. 

The plan of the thesis is as follows. In 
Chapter II, a brief review of the literature on heuriS' 
tic procedures for flow-shop scheduling problem is 
presented. Chapters III is devoted to the development 
of basic frame-work which is used for structuring the 
proposed heuristic algorithm presented in Chapter IV. 
The proposed algorithm has been tested on standard 
problems available in the literal ure and also on a 
large number of randomly generated problems of varied 
sizes. The results are presented in Chapter V. 



CHAPTER II 


LITERATURE SURVEY 

In the area of flow-shop scheduling, considerable 
research work has been reported in the literature. 
Minimization of make-span time and average flow time are 
generally considered as the optimization criteria. A 
thorough survey of the literature indicates that the 
flow shop problems have been tackled using mathematical 
programming, complete enumeration, Monte Carlo sampling, 
decomposition and Heuristic approaches. Since the pre- 
sent work is primarily concerned with the development of 
heuristic approaches for flow-shop scheduling problems, 
only the literature relevant to this is reviewed in the 
following section. 

2.1 Heuristic Approaches ; 

Heuristic models are developed utilizing logic 
and common sense derived from the observation of the 
behavior of solutions under varied conditions. Heuristic 
procedures never guarantee optimality. These are some 
sort of intuitive proposals that would give, in most of 
the cases, near optimal results. The effectiveness of 
a heuristic is judged by its easiness and closeness 
of its results to the optimality. 
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A number of heuristic models have been proposed 
by different researchers. The notable among them are 
due to Palmer (10), Gupta (3), Petrov (ll), Cambell, 

Dudek and Smith (2), Hobeshima (7), Krone and 
Steiglitz (6), Gupta and Maykut (4). Palmer has sugges- 
ted a quick method of scheduling. He assigns priority 
to items by determining a numerical slope indoz for each 
job and constructing a sequence in non-decreasing order 
of the slope indez. 

In Gupta's method (3), a numerical value is 
assigned to each job on the basis of a predefined func- 
tion and the sequence is constructed in the ascending 
order of the values. Both Palmer and Gupta have used 
the analogy between scheduling and sorting problems, but 
Gupta claims that his method is superior to Palmer’s 
algorithm. 

Cambell, Dudek and Smith (2) have treated the 
problem by constructing (m-l) sub problems, where m is 
the number of machines. For each sub-problem, a sequence 
is obtained using Johnson's 2-machine optimum method. 

The sequence which jrields the minimum make-span is 
finally selected. 

Petrov's algorithm (ll) is similar to Campbell, 
Dudek and Smith's algorithm (2) with one important 
difference that it does not use Johnson's algorithm. 



8 


Petrov's algorithm has been used in the present work 
for the generation of initial sequence. The details 
of this algorithm are presented in Appendix A. 

Krone and Steiglitz (6) have proposed a heuristic 
programming approach with mean job completion time 
objective function. They fix a starting schedule by a 
pseudo-random algorithm and use a heuristic programme 
which computes a number of trial sequences. The 
algorithm terminates when local optimality of the 
current trial solution is verified by searching its 
neighbourhood exhaustively. 

Page (8,9) has developed four heuristics based 
on sorting techniques, merging (M), pairing (P), indivi- 
dual exchanging (IE), and group exchanging (GE). The 
merging approach considers the list of jobs as a collec- 
tion of strings of items. Initially the list is divided 
into n strings, each containing one job. Each successive 
pair of strings is then merged into a single ordered 
string containing two jobs. The ordering is based on 
the minimum make-span for the pair of jobs. Continuing 
in this manner, the number of jobs per string is increased 
while the number of strings is decreased until a single 
ordered string of n items remains. The pairing proce- 
dure is almost identical to that for merging, with the 
exception that once a string has been formed, adjacent 
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jobs will remain permanently adjacent in any larger 
strings developed. 

The individual exchanging heuristic is based 
on the exchange method of sorting. Starting with a 
given job order, a test is performed to see whether 
each successive pair of adjacent jobs should remain 
as they are ordered or exchange positions. An exchange 
is made whenever the new sequence has a lower make-span. 
If an advantageous exchange is found on the first pass, 
a second pass is ma.de. The process continues until a 
pass fails to yield a profitable exchange. 

This is infact a very inefficient way of making 
the exchanges. Some of the exchanges can be dropped 
outrightly because they would worsen the sequence if 
made. Moreover, Page's procedure confines only to 
adjacent pairwise exchange and the higher order exchanges 
are not considered at all. In the present work an algo- 
rithm has been developed which takes care of these draw- 
backs. 

lastly, Page’s Group Exchange Method is quite 
similar to individual exchange method except that groups 
of jobs are exchanged rather than individual jobs. 



CHAPTER III 


FOUNDATIONS OE PROPOSED METHODOLOGY 

In. Chapter II a review of the heuristic procedures, 
available in the literature for flowshop scheduling 
problems, was presented. It was pointed out that, in 
general, these heuristic procedures are computationally 
efficient but do not necessarily yield optimal sequences. 
The problem under consideration is to evolve a methodology 
which when used on an arbitrary sequence would improve 
upon the sequence which may further be used as an input 
sequence for further improvement. The procedure would 
have to be iterative in nature and would terminate when 
no further improvements are feasible considering a pres- 
pecified measure of effectiveness. In the present work 
makespan time has been used as the criterion of effective- 
ness of a sequence. The initial sequence, instead of 
being generated arbitrarily, ia developed using one of the 
best heuristics available in tbe literature for flowshop 
scheduling problem. 

The problem can now be formally stated as follows; 
3.1 Statement of the Problem ; 

Given a sequence for a static flowshop scheduling 
problem which has been generated arbitrarily 
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or through, the use of a heuristic, develop a methodology 
which exploits certain characteristics of this schedule 
for the generation of improved sequence considering 
makespan time as the criterion of effectiveness.' 


The above problem statement actually focuses on 
the ultimate goal, the statement would acquire the 
following structure if it is analysed in terms of how 
to achieve the goal. 

Consider a static flowshop involving n jobs and 
m machines. Let t^ represents the processing time of 
j-th job on i-th machine. Further, let 

S = J(l)-J(2) .... - J(j) J(n) 


represents a sequence with J(j) representing the job 

in position j of the sequence. Given a sequence S and 

processing time matrix T = ft^j , find the job 
t ft ^ 

positions j- ^ and jV^ for the cyclic pairwise exchange 
in the given sequence such that the makespan time of 
the new sequence S*", obtained after making the said 
exchange, is lower than that of previous sequence. The 
term cyclic pairwise exchange is explained in the next 
section. 


Each exchange constitutes an iteration, if an 
iteration yields an improvement, the next iteration is 
undertaken. The process is continued until no further 
improvement is possible. 



3.2 Cyclic Pairwise Exchange ; 

Consider a sequence S = J(l)~J(2)... J(n) 
for cyclic pairwise exchange which is to be effected 
to displace J(k) job to position i. The exchange under 
consideration can be specified as fj(k), i'j . The pro- 
cedure involves the removing of job J(k) and shifting 
the string of remaining jobs lying between positions 
i and k by one position towards k such that the i-th 
position is vacated and the k-th position gets occupied. 
In this process the job J(k) gets placed at position i. 

It should be noted that jobs beyond i-k range remain 

unaffected. 

Por the sake of illustration, consider a 
sequence S = 4-3-1-2-5. let us say job J(4) = 2 is to 

be displaced to position i = 2, in other words the 

exchange j 2,2 j is to be made. The new sequence after 
making the said exchange will be ; 

S* = 4-2-3-1-5, 

3 • 3 An Approach for the Development of Proposed 
Methodology ; 

In this section the basic framework which is 

of 

used for the development^ the proposed methodology is 

jt" + is the job position in the sequence from where 
0 0 the job should be taken out. 
t f 

jr is the job position in the sequence where the 
job J ( J ou -£ ) has to come in. 
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presented. Throughout this work an operation is speci- 
fied by task position in the sequence and the machine on 
which it is to he performed. For example, an operation 
represented as (i, j) connotes that the job in the j-th 
position of the sequence is to he performed on machine i. 

As has been pointed out earlier, certain important 
characteristics of the given sequence which is assumed 
to be sub-optimal in character, need to be exploited for 
the identification of cyclic pairwise exchanges. Obviously, 
one such charact eristic would he the idle time associated 
with some of the jobs and machines. In the process, a job 
may have to wait for want of availability of the machine 
and similarily a machine may have to wait for want of 
job. These waiting times arc referred as idle times 
here. For the identification of such idle times, one 
alternative is to plot Gantt chart for the given sequence. 
However, to facilitate the computerization of the proposed 
methodology a matrix approach is adopted for obtaining 
the idle times. A matrix referred as Earliest Finish 
Gantt (EFG) matrix is developed first and is used for 
generating the idle times for the various operations in 
a matrix form. 
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3.3.1 earliest Finish Gantt (BEG-) Matrix ; 

This is an m x n matrix which represents the 

earliest completion times of various operations and is 

denoted by G- = /g. .1 . The element g. . indicates 
J l & ij/ . 13 

the earliest completion time of operation (i,j). 

The following relationships can be used for 
the development of EFG matrix; 


S, 


'ij 

= max 

L(g ij 3-1* 

S i--1.3‘ 

'0, 3 

= 0 

3 = 1, 

p * « # n 

o 

os 

•H 

= 0 

i = 1, 

2 rn 


It should be noted that the element g^ n will 
represent the makespan time of the given sequence. The 
EPG- matrix information can be presented in the form of a 
chart which will be referred as EFG chart. On this 
chart, corresponding to each operation there is a hori- 
zontal bar which is called Operation Block. 


3.3.2 Idle Time Matrix ; 

The idle time matrix 1 = 

by the following relationships; 



can be generated 


I i, 3 S i-1, 3 “ g i ? 3"1 

I*l * * =: 0 t 2. 2 » O m f Yl 

-L # J 

I i . i = _i . i i = 2,3,..., m 



An element of this mafias can take any value. 

Negative I. • would imply that the job in the j-th 

position of the sequence S awaits the i-th machine 

for I. . time. A positive I. . signifies that the 
l J J t ? J 

i-th machine waits for the jot 1 in the j-th position for 

this period. Obviously, I. . = 0 means that the job 

1-9 J 

in position j is immediately taken up by i-th machine. 


The concepts of latest Finish Gantt (LFG) matrix, 
slack time (matrix) and critical path have been intro- 
duced to facilitate the development of the methodology. 

3.3.3 latest Finish Gant t ( IPG ) Mat rix ; 


This matrix represents the latest completion times 
of the various operations and is denoted by G' = j* g^ .. j , 
where g'. . indicates the latest completion time of 

3- 9 J 

the operation (i, j). It may be noted that for a sequence, 
g! . will always be greater than or equal to g. .. An 

1 9 3 J 

algorithm has been developed for the generation of LFG 
matrix. 

The various steps of the algorithm are given 

below ; 

Step 1 ; For a given sequence S = J(l)-J( 2) , . . . , J(n) and 
machine ordering l,‘2,... s m, obtain the EFG matrix 
G = / g^^ | . Find out the makespan time M g = g^ n « 
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Step 2 ; Reverse the sequence as veil as the machine 
ordering, new sequence will be 8 p = J(n)-J(n-1) „ . . . , J(l) 

and the machine ordering m 5 m-1, . .., 2,1. 


Step 3 ; Rind EPG matrix for the new case. Let this 


matrix be G- = f 

r ~ 


where g. . = max <* g. . , g. + t . n 

& i3 j h, 3“1 & i-l s 3 I m+l-x,n+l -3 


g n = 0 

o » j 


3 ■ — 1j^s'*'j n 


g. 


i , 0 


= 0 


r — _L 4 2 a ...c m. 


4s Obtain a matrix G = 


( s il 


where , 


13 


KL 


5 i3 


% Reverse the matrix G in such a manner that 

i-th row becomes .(m+l-i) th row and j-th column becomes 

= ? 

(n+l- 3 ‘)-th column. Let this matrix be called G . The 
G' matrix gives the latest start time of various opera- 
tions . 

Step 6 s Obtain the LPG matrix G’ = J gl „ ; 1 , 

— ^ x, j J 


where 


"ij s ij + ^,1(3) 


3.3.4 Slack Time ; 

Slack time is the maximum amount by which an 
operation can be delayed without affecting the makespan 
time. The slack time matrix K = f can be obtained 
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by subtracting the EFG- matrix from the LFG- matrix. Thus 
K = G* - G 

It needs to be pointed out that the operations 
(m, n) and (l, 1) will always have zero slack time. 


As would become obvious in Chapter IV, the slack 
time matrix becomes necessary for determining the j. 
position in the given sequence for making the cyclic 
exchange. Instead of generating the slack matrix for 


the complete sequence, only the matrix for the sequence 
of size (n-l), obtained after removing job J(j Q ^) from 
the sequence, need to be generated. 


3 , 3.5 Critical Path ; 

Critical path is the path obtained by joining 
those (m+n-1) operations each of which has zero slack 
time. Obviously, the makespan time is changed if any 
of the operations on the critical path is shifted. There 
can be more than one critical path corresponding to a 
given sequence. 

A critical path is denoted by a string of (m+n-1 ) 
operations. For example, (3,4)-(2,£)-(2,3)-(2,2)-(l,2)-(l,l) 
can be a critical path for a 3-machine, 4- job problem. 

The critical path(s) corresponding to a given 
sequence can conveniently be depicted on a grid diagram. 

In a grid diagram, the horizontal and vertical lines 



denote the machines and jobs, respectively. The nodes 
represent the operations. In Fig. 3.1 a grid diagram 
for a 4-machine, 6-job problem is presented. Further, 

JL_ 

2 

3 

machines ~ ' 

4 


Job positions 

Fig. 3.1 

the critical paths 

(4,6)-(4,5)“(b,5)“(3,4)“(2,4)-(2 f 3)-(2,2)-(l,2)-(l,l) 

and (4,6)-(4>5)-(4,4)-(4,3)“(3»3)-(2,3)-(2,2)-(l,2)-(l,l) 
are indicated on this diagram. 

Following are some of the important properties of 
critical paths ° 

1. If an operation has zero slack time, at least one 
critical path will pass through it. 

2. All critical paths start from operation (m,n) 
and end at operation (l,l). 

3. If some operation (i,j) belongs to a critical path 
then the next operation on that path will either be (i-l,j) 
or (i, j-l) provided (i,j) is not (l,l). 
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4. There will not. be any gaps between the operation 
blocks lying on a critical path. 

5. If the slack time of operations (i-l,j) and 
(i,j-l) is zero, the idle time associated with the opera- 
tion (i,j) will be zero. 

6. If the operations (i-l,j) and (i,j) have zero slack 
then the operation ( i » j ) will have non-negative idle 
time. 

7. If the operations (i, j-l) and (i,j) have zero 
slack then the operation (i, j) will have non-positive 
idle time. 


Two important theorems along with their proofs 
are presented next. 


Theorem I s If an operation (i,j) has zero slack then 
at least one of the two operations, viz., (i-l,j) and 
(i, 3 ~l) will have zero slack. 

Proof ; In the SPG- matrix, the earliest finish time of 
an operation (i,j) is given by 






max 




1 s 3 




3-1 j 


i + t_ 


J(j) 


The earliest start time for an operation (i,j) can be 
obtained bjr subtracting t. Tf • \ from g. .. Thus the 
earliest start time is given by, 





g i, 3-1 



max 



The above relationship suggests that the earliest 
finish time of at least one of the two operations, viz. , 
(i-1, j ) and (i, j-l) would coincide with the earliest 
start time of operation (i,j). The operation whose ear- 
liest finish time coincides with the earliest start time 
of operation (i,j) can not be shifted relative to 
operation (i,j) which already has zero movability because 
of its zero slack. Therefore, at least one of the 
operations (i--l, j ) and ( i 5 j -1 ) would have zero slack. 

Theorem 2j If an operation (i,j) has zero slack then 
at least one of the two operations, viz., (i+l,j) and 
(i,;j+l) will have zero slack. 

Proof ; The movability of operation (i, j) is dictated by 
the operation (i+l,j) and (i,j+l) which are the only 
direct successors to operation (i,j). By contradiction, 
if both the operations (i+l,j) and (i,j+l) have positive 
slack then the operation (i,i) will have no operation 
restraining its movement towa,rds the greater value on 
the time amis. Since the operation ( i , j ) has zero slack, 
it implies that at least one of the operations (i+l,j) 
and (i,j+l) will have zero slack. 

The above theorems essentially imply that it is 
always possible to construct at least one critical path 
for a given sequence. 



3.3. 5.1 Algorithm for Critical Path Generations 

An algorithm has been developed for generation of 
critical paths corresponding to a given sequence. In 
order to save on computational effort the algorithm uses 
the information about Idle time matrix rather than 
slack time matrix. The various steps for the algorithm 
ares 

Step I s Consider the Idle time matrix I = / 1^ 
for the given sequence. 

Step 2s Start from the operation (m,n). Set i = m, j = n 
and 0=0, where C identifies the critical path number. 

Step 3 s Set C = C + 1 

Step 4; If I 4 • is negative, go to Step 5 

3- 9 J 

positive, go to Step 6 . 
zero , go to Step 7 

Step 3 ° Connect operations (i,j) and (i,j-l). Set j = j— 1 
Go to Step 8. 

Step 6 s Connect operations (i,j) and (i~l,j). Set i = i-1 
Go to Step 8. 

St ep 7 s Connect operation (i,j) with (i,j-l) and ( i-1, j ) 
separately. Set, i = i-1. Keep |(i, j-1),* C \ in the 


memory. 
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If i = 1 or j = 1, go to next Step . Otherwise go 


to Step 4. 


Step 9; If i = 1, join (1, j) with (1,1) decreasing j 
by one in each step. 

If j=l, join (i,l) with (1,1) decreasing i 
by one in each step. 


Step 10 ; If there is nothing left in the memory, STOP. 

Otherwise, take the first set <j ! (i,j)j cl from the 

t J 

memory and remove it :from there. Join operation (i,j+l) 

* 

with (m,n) through critical path C, and go to Step 3. 


3.4 


The various concepts introduced in this chapter 
are illustrated through an example. 


Consider a 4-machine, 6- job flow shop with the 
following processing time information; 

5 10 6 6 3 9 

16 3 9 4 5 7 

T = 

12 4 8 7 11 6 

3 7 5 13 3 8 


In all, there will be 61 feasible sequences, 
let us consider one of these feasible sequences. Let 
this sequence be S = 4-5-1-6-3-2. In this sequence 
J(4) = 6 indicates that job 6 takes 4th position in the 
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sequence and will be processed when the first 3 jobs, 
viz., 4-5-1 have been completed. 

EP& Matrix and Char t; 

for the above sequence and the machine ordering 
1,2, 3,4 , the EPG matrix will be 


6 

14 

19 

28 

34 

44 

10 

19 

35 

42 

51 

54 

17 

30 

47 

53 

61 

65 

30 

33 

50 

61 

66 

73 


The corresponding EPG— chart is shown in Pig. 3.2. The 
shaded areas indicate the machine waiting times. 



The idle time matrix will be as follows; 


0 

0 

0 

0 

0 

0 

6 

4 

0 

-10 

-11 

-10 

10 

2 

5 

- 5 

- 2 

- 7 

17 

0 

14 

3 

0 

- 1 
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Here, ^2 4 = means that the operation (2,4) 
awaits machine 2 for 10 units. 

Note that the various shaded areas in Nig. 3.2 
correspond to the positive numbers in the I matrix. 

IPG Matrix : 

The algorithm given in Section 3.3.3 is applied 
for generating the LPG matrix. 

Step I s Prom the G-matrix, already obtained, M g = 73. 
Step 2 ; S^ = 2-3-6-1-5-4 , and machine ordering: 4,3,2,! 


G = 


G 


Step 5 ° 


G’ 


7 

12 

20 

23 

26 

39 

11 

20 

26 

38 

49 

56 

14 

29 

36 

54 

59 

63 

24 

35 

45 

59 

67 

73 _ 

66 

61 

53 

50 

47 

34 

62 

53 

47 

35 

24 

17 

59 

44 

37 

19 

14 

10 

49 

38 

28 

14 

6 

0 

0 

6 

14 

28 

38 

49 

10 

14 

19 

37 

44 

59 

17 

24 

35 

47 

53 

62 

34 

47 

50 

53 

61 

66 
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6 

14 

19 

37 

44 

59 

14 

19 

35 

44 

53 

62 

24 

35 

47 

53 

61 

66 

47 

50 

53 

61 

66 

73 


Slack Time Matrix ; 


Using the matrices G* and G, the slack time matrix 
K is obtained, 

0 0 0 9 10 15 

4 0 0 2 2 8 

t r 

Xi. = 

7 5 0 0 0 1 

17 17 3 0 0 0 

Criti c al Path ; 

The various critical paths are constructed using 
the algorithm given in Section 3.3.5. The basic input 
is the I matrix which is reproduced below; 

0 0 0 0 0 0 

6 4 0 -10 -11 -10 

10 2 5 - 5- 2- 7 

17 0 14 3 0 -1 

Critical path C = 1; 

Starting from the operation (4,6) and using the 
I-matrix we get, 
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4,6 

C 4,5 


I* 


3,5 

3,4 

3.3 

2.3 


- 1 , so move horizontally to ( 4 , 5 ) 

0, move to ( 3 , 5 ) and keep 1 ( 4 , 4 ),- 1 \ 

* j 

in the memory 
-2 ? move to ( 3 , 4 ) 

- 5 , move to ( 3 , 3 ) 

+5, move to (2,3) 

0, move to ( 1 , 3 ) and keep j ( 2 , 2 ); 1 j 


in the memory 

Since i = 1, join (1,3) with (1,1)- 

Hence the critical path C = 1 is: 

(4»6)-(4,5)-(3,5)-(3,4)-(3,3)-(2, 3)-(l,3)-(l,2)-(l,l). 

look into the memory : We have f(4,4),- 1*1 and j(2,2)$ lj. 
let us take the first set first for the generation of 
another critical path and remove it from the memory. 

Critical Path C = 2i 

Take (4,6) to (4,5) portion of critical path C = 1. 
Start from ( 4 , 4 ), 


II 

H 

+ 3, 

move 

to (3,4) 


I 3,4 

- 5, 

move 

to (3,3) 


I 3,3 = 

+ 5, 

move 

to (2,3) 

/( 2 , 2 ), 2 ) 

II 

C\J 

H 

0 , move to 

( 1 , 3 ) and keep 

memory. 

Since 

I = ! 

1, join(l,3) with (1,1). This 


completes the identification of the second critical path- 
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Thus, the critical path C = 2 is ; 

(4.6) -(4,5)-(4,4)-(3,4)-(3,3)-(2,3)-(l,3)-(l,2)-(l,l) 

how the memory has ^(2,2),* lj and|(2,2)j 2j sets., 

Selecting |( 2,2), 1] first, we generate next critical 
path. 

Critical Path C = 3s 

Portion of critical path C = 1 from (4,6) to (2,3) is; 

(4.6) -(4,5)-(3,5)-(3,4)-(3,3)-(2,3) 

Starting from (2,2). 

Ij 2 = + 4 , move to 1^ 2 

Since i = 1, join (1,2) with (1,1) 

This completes the identification of 3rd path, viz., 

(4.6) - i (4,5)—(3,5)-(3,4)-(3»3)~(2,3)-(2,2)-(l,2)-(l,l) 
ho more addition in the memorjr. 

Taking the only set |( 2,2), 2^ left in the memory. 
Critical Path C = 4; 

Portion of critical path C = 2 from (4,6) to (2,3) is; 

(4.6) -(4»5)-(4 # 4)-(3,4)-(3,3)-(2,3) 

Start from (2,2) 



1,2 


+ 4, move to I 



Since i = 1, join (1,2) with. 
Hence the critical path 0=4 
(4,6)~(4,5)-(4,4)-(3,4)-(3»3) 


Nothing left in the memory 


( 1 , 1 ). 

is 

-(2,3 )“•(!» 2)-(l,l) 


Therefore, STOP. 



CHAPTER IV 


METEOROLOGY POP, IMPROVING THE SEQUENCE 


4.1 Introduction ; 

In this chapter a method ology for improving 

a given sequence considering makespan time as the 

criterion of optimization is presented. The methodology- 

draws upon the conceptual framework of slack matrix and 

critical path which were introduced in Chapter III. Por 

improving upon a given sequence, the methodology uses 

cyclic pair wise exchanges for which j ou ^ and 3^ need 

to he identified. The identification of i. and 3 . 

un °out 

has to he such that the exchange results in reduction 
of the makespan with respect to Hie initial sequence. 

The process of cyclic exchanges is carried out till no 
further improvement is possible. The methodologies for 
the identification of 3 , and 3. are discussed separa- 

tely and then the complete algorithm is presented. 


4.2 Identification of 3 ^ Position ; 

Por the identification of 3 ^ position the 

following heuristic rule is developed. 


Rule; Only those jobs should he considered for removal 
for the cyclic exchange whose positions appear more than 
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once in at least one of the critical paths corresponding 
to the given sequence. 

4.3 Identification of .ih Po s ition s 

A systematic methodology has been developed for 
the purposes of identifying the j. position. The metho- 
dology utilizes the concept of 'unsatisfied* critical 
path which is discussed and explained with the help of 
an example. 

4 . 3.1 Unsatisfied Critical Path ; 

A critical path is said to be unsatisfied by an 
- 

exchange |J(k), ij if; 

1 . The 30b positions i and k appear 3'ust once in cri- 
tical path under consideration, and 

2 . The job positions i and k have same machine 
number in their respective operations of the critical 
path under consideration. 

If any of the above stated considerations is 
violated, the critical path is said to have been satis- 
fied by the exchange. 

To illustrate the above definition an example is 
presented . 

Let us consider a critical path which is represented 
as; (4 ? 6)-(4,5)-(3,5)-(3,4),(3,3)-(3,2)-(3 ,D-(2,l)-(l,l). 



For this critical path it is observed that the j oh 
positions 2 and 4 appear only once and for both the 
job positions the associated machine is the same, i.e., 

machine number 3. Therefore the cyclic exchanges 

f f ^ 

jj(4),2j and jj(2), 4j will ' unsat isfy' the given 

r ^ 

critical path. On the other hand the |J(3), 5? is one 
of the exchanges that would satisfy this critical path. 

Having introduced the concept of 'unsatisfied' 
critical path, a theorem is presented for the identifi- 
cation of plausible positions. 

Theorem ; If a cyclic exchange j J (k), ij results in 
some of the critical path(s) of a given sequence to 
remain ’unsatisfied', then the new sequence obtained 
by the said exchange will have makespan greater than 
or equal to that of the original sequence. 

Proof; Since the exchange fj(k), il unsatisfies some 

of the critical paths, both the operations having job 

positions i and k and lying on these critical paths, 

will have the same machine number. Since there are no 

gaps between operations lying on a critical path, the 

operations lying between job positions i and k on the 

critical paths can only have non-positive idle times. 

r 

Therefore, if the exchange <J(k), i is made, the time 
span contribution due to the operations lying between 
the two operations which took part in the exchange will 

l I T. KANPUR 

CEN I « At. ? f 

Now * * 



not decrease in any case? hence the makes pan will not 
decrease. However, the above exchange may cause the 
other operations to arrange themselves in such a manner 
that some gaps get introduced on the critical path(s) 
under consideration. Such an exchange will increase 
the make span time. 

The ah ove stated theorem is explained below with 
the help of an example. 

Consider a critical path (4»6)—(3»6)-(3»5)-(2,4)- 
( 3 , 3 )-( 3 » 2 )-( 2 , 2 )-(l, 2 )-(l,l) which is depicted in a grid 
diagram as well as in EPG chart as shown in mg. 4.1 


J ohs 



Pig. 4.1 


|j( 3),5j is one of the cyclic exchanges that would 
unsatisfy the above mentioned critical path. If this 
exchange is made, the time span between operations (3,3) 
and (3,5) will not decrease, because there are no gaps 
between the two operations which could be reduced or 
vanished completely. However, some gaps can get 
introduced on the critical path due to reshuffling of 
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operation blocks. Hence the makespan will either 


increase or remain same. 


The theorem discussed above gives rise to the 
following rule. 

Rule ; For improving a given sequence, only those 
cyclic pairwise exchanges should be considered which 
satisfy all the critical paths corresponding to the 
given sequence. 

It should be noted that corresponding to an 
identified j QU+ there can be number of candidates. 
The procedure to determine the set of jobs which can 
be removed^ and their corresponding sets of j. candi- 
dates is explained with the help of an example below. 


net us say following are the only critical paths 
corresponding to a given sequence 

S = J(1)-J(2)-J(3)~J(4)-J(5)-J(6) 


Critical 
path 
number C 


Critical Path 


(4.6) -(4,5)-(3,5)-(3.,4)-(3,3)-(2,3)-(l,3)-(l,2)-(l f l 

(4.6) -(4,5)-(4,4)-(3»4)-(3»3)-(2,3)-(l,3)-(l,2)-(l,l 
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Apply the rule mentioned in Section 4.2 to deter- 
mine the set of jobs which can be removed from the given 
sequence. We get , 

J(5) because job position 5 appears twice in critical path 3. 

J(4) because job position 4 appears twice in critical path 2 

J(3) because job position 3 appears thrice in both critical paths 

hext step is to determine the sets of j in candidates 

for each of J(5), J(4) and J(3). 

Consider each of these three j obs separately, say 
job J(5) is considered first, how apply the rule mentioned 
in Section 4.3 to determine the two sets of j^ n candidates 
with respect of the two given critical paths. We get, 

Critical path number C Set of candidate s 

1 all positions 

2 1 , 2 , 3 , 4 

The intersection of these two sets, i.e. , 1,2, 3, 4 

is the actual set of j^ n candidates corresponding to J(5). 

Hence, J(5) can be placed in position 1, 2,3,4 of the 
given sequence S. 

In the similar way, the sets of j^ n candidates 
corresponding to jobs J(4) and <T(3) can also be deter- 
mined. 

Once the set of j. candidates is identified, the 

111 

problem remains to select a specific j. for exchange with 
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the already determined i 

'-’our 

4 .3. 2 Selection of 


Once 3 ' , is specified tlie information from slack 

and idle time matrices is utilized for the selection of 

a specific j in out of the set of its candidates. The 

procedure involves the selection of that 3 . which results 

m 

in maximum reduction in the makespan time. This implies 
that one would have to generate the EFG- matrix for 
calculating the makespan time after each exchange is 
made. This can he fairly time consuming. Therefore, 
procedure has been developed which determines the change 
in the makespan time due to an exchange { -J( 3 0U ^. ) •, j^ n '1 
without generating the EFG matrix. The development of 
the procedure is explained in the following section. 


4. 3.2.1 Procedure : 

Given a sequence S = J(l) -1(2) ... J(n) 
let 3 in and 3 0 ut corrGS P° nG 'fc 0 3 and 3 ' s respectively. 
After removal, of job J(j ! )9 the sequence will have (n-l) 
jobs. Corresponding to this sequence of (n-l) jobs EFG , 
idle time and slack matrices are generated. Fete that 
these matrices will be of size m x n-l. 

If the removed job J(j ! ) is placed at position j 
in the sequence of (n-l) jobs then what would be the 
increase in the makespan, is determined using a proce- 
dure discussed below. 
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At this stage we have the sequence of (n-l) jobs 
for which an EFG- chart can be developed and further the 
makespan can be evaluated, let us say that the 
removed job is to be placed in position j. This would 
displace the present job in position j along with its 
followers by one position. In Fig. 4.2 a portion of the 
EFG- chart which corresponds to two jobs, viz., J(j-l) 
and J(j) is presented. The positive idle times are 
shown by hatched lines. One can utilize these idle times 
advantageously for accommodating the job to be inserted. 

_ JM; _ 2 

M/c ! g i,j-l j g l,j 

M/c 2 ^2, j~l §2, j ! 



M/c (m-2) 
M/c (m-1) 
M/c m 



Fig. 4.2 ; EFG Chart 

Further, the slack times associated with job in position 
j can also be utilized ..for this purpose. Therefore, the 
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total available gap at machine i for the insertion of the job 
J(j’) is given by. 


+ thj- 01 

However, there is a possibility that the actual 
available gap is some what less. The reduction in the 
gap will occur when the operation block of job J(j') 
corresponding to machine (i-1) is large enough to 
cover the operation block (i, j-l) as well as a portion 
of the total gap mentioned above. This kind of situation 
is depicted in Pig. 4.3. 


M/c- i-1 
K /c i 


; T :i-1, J(j' ) 


Si, j-i 


Pig. 4.3 

Thus the net gap available for accommodating 
the operation of job J(j’) on machine i, i.o., t^ j( j 1 ) i- s 

P i = k i s j + max f X i,j s0 ^ “ maz C( T i“Si 3 -«.i)»°] 

where, x ± == + t^-^ ) and ^ = 0 

If t. rf -\ is greater than P - , the available gap would 
have to be increased. let the minimum shift which 
needs to be given to the jobs lying after position (j-l) 
to create this additional gap required corresponding to 
machine i be x^. This implies that 
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x ± = *±,xy) - k i,j " max [I ±,j’ 0] 

+ max [(r ± - g ± 0] 

So far we have been confining our attention towards 

accommodating i-th operation of J(j ? ) on machine i only. 

Since this job is processed on all the machines, the 

minimum increment required in the makes pan, M*, will 

be Max (x^). 
i 

4 . 4 Alg orithm ; 

The following are the steps of the final algorithm 
developed for improving upon a sequence. 

Step Ij Obtain a sequence S = J(l)-J(2) . . . J(n) using 
any one of the known heuristics. 

Step 2; G-enerate EDG-matrix G = (g. . | , idle time matrix 

“ • * . A * 1 1 J ) 

I = J I. . l < for the sequence S. Determine makespan M c< . 
t i > 3 I ° 

Step 3 i Construct critical path(s) using I matrix. 

Step 4 : Store the job positions which appear more than 
once in' any of the critical paths that have been cons™ 
tructed in Step 3. let these job positions are stored as 


set 


( V(k) 


k = 1,2 


... X j. 


K is total number of such 


positions. 


Set k = 1 and M g =oC 


Step 6 i Set v = 7(k) 



Step 7; Obtain, a sequence S' = J’ (l)-J' ( 2 ) . . . J* (n-1) after^ 
removing J(v) from sequence S. 

St ep 8 s Generate EPG-matrix , I -matrix, 2-matrix with 

I 

respect to sequence S'. Determine make span M 0 . 

s 

Step. ,9..; Identify the set of • candidates - in sequence S 
corresponding to job J(v). Let this set be denoted as 

; "7 

1 H(q); q. = 1,2,... L j , where L is total number of 
candidates corresponding to J(v). 


Step. 10: Determine increment in make span for each candidat 

H(q), Pick up the candidate which incurs the minimum 

increment and store it in i. . Let the minimum increment 

111 

be called HIM. 

Ste p 11 : If (HIM + M ) is less than M , go to Step 12. 
Otherwise go to Step 13. 

Step 12 : Set IlCfllT = HIE and j Q t = v. Determine 
m" = Ii' + IMCMT . 

b B 


Step 13: If k = E, go to Step 15. Otherwise go to Step 14. 
St ep 14; Set k = k+1 and go to Step 6. 


Step. 15 g If INCH! is greater than or equal (M g 
to Step IQ. Otherwise go to Step 16. 

Step 1 6; Make cyclic exchange | J( J out ), J iri j 



on sequence S 


, , ■ « 

to Occam a new sequence g . 

Step. 17: Set S = 8*' and go to Step 2. 




! Apply Petrov T s Alg* to obtain an Initial S 
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Step 10 ; Print , *Fo more improvement possible on the 
original sequence 1 , Hence STOP. 

A general logic flow diagram of the above stated 
algorithm is presented in Pig. 4.4. 

4 , 5 Ex ample ; 


To illustrate the various steps of the algorithm 
presented in the last section an example follows; 

Given a processing time matrix 


10 

12 

1 

13 

6 

5 

19 

5 

7 

15 

7 

8 

19 

5 

14 

16 

19 

18 

7 

10 

11 

16 

17 

11 

7 

12 

13 

21 

8 

12 

3 


13 

9 

11 

15 

18 

3 

14 

14 

11 

3 

11 

15 

13 


for a 5 ™ job, 3 -machine flow shop problem. 

Step 1; Using the Petrov’s algorithm the following sequence 
is generated. 


8 = 

3 - 

5 - 

1 - 

8 - 

- 9 - 

.. 4 . 

- 7 - 

- 6 - 

- 2 


1 

7 

17 

22 

29 

42 

61 

6 6 

78 


g 

14 

32 

51 

69 

88 

104 

118 

125 

Cf = 

20 

37 

44 

63 

82 

*sj“ 

o 

1 — 1 

111 

129 

139 


32 

38 

65 

76 

97 

107 

120 

142 

150 


46 

57 

83 

98 

111 

125 

136 

145 

153 
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I 


0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

-2 

3 

.-10 ■ 

-22 

-27 

-27 

-38 

-40 

9 

-6 

-5 

7 

6 

6 

0 

7 

- 4 

20 

5 

6 

- 2 

6 

7 

4 

9 

- 3 

32 

O 

— O 

o 

O 

- 7 

- 1 

- 4 

- 5 

6 

5 


Make span M g 


153 


Step 3. ° There will be onlj? - one critical path as given 


(5,9)-(4,9)-(4,8)-(3,8)-(2,8)-(2,7)-(2,6)-(2,5)- 

(2,4)-(2,3)-(l,3)-(l,2)-(l,l) 


Step 4 ; The ■ 

total 

number of 

can did at 

es for 

^ out 

is 

and they are 

7(1) 

= 3 

, 7(2) 

= 8 

and 

7(3) = 

: 9. 


Step 5 ; k = 1 

and ? 

_U 

0 c 






Step 6 : v = V(l) = 

3 







Step 7* Remove J(3) 

1 = 

1 , we 

get 





II 

m 

- 5 - 

- s 

- 9 - 

4 - 

7-6 

- 2 



Step_ 3 :: 

r i 

7 

12 

19 

32 

51 

56 • 

68 


i 

9 

14 

33 

51 

70 

86 

100 

107 

G = 

20 

37 

49 

64 

06 

93 

111 

121 


32 

38 

60 

79 

39 

102 

124 

132 


46 

57 

75 

92 

106 

117 

127 

135 


p 

0 

0 

0 

0 

0 

0 

0 

0 


1 

-2 

-2 

-14 

-19 

-19 

-30 

-32 

I = 

. 9 

-6 

-4 

2 

6 

0 

7 

- 4 


20 

5 

11 

4 

7 

4 

9 

- 3 


32 

-3 

3 

4 

-3 

-4 

7 

5 


slow 


3 



-cj 






0 

0 

a 

2 

14 

19 

19. 

30 

39 






0 

0 

Q 

0 

0 

0 

0 

7 




Jf 

= 

14 

14 

14 

12 

7 

7 

0 

5 






! IS 

27 

16 

12 

13 

9 

0 

0 





t 

19 

IS 

16 

12 

12 

12 

5 

0 

1 



Ms 

Ice span H,!. 

= 

135 







Step 

.9” 


All 

CD 

v 0 

3°; 

0 posit. 

ions a 

ire 

the 1 

c and idat 

es for 

joh 

J(3 

) 

= 1. 










Step, 

JLO 

• 

Inc 1 

cement 

s correspo: 

a ding 

to 

each 

j . Candida 
J in 

.te 

are, 















Joh position: 

1 2 

3 

4 

5 

6 7 

8 

9 



IncroT! 

I — * 

CD 

P 

CP" 


24 17 

13 

15 

15 

15 15 

21 

36 


Candidates 4 , 5 , 65 ? yield minimum increment. So any 


of these should he selected, 
selected. 

Therefore, jh r = 4 and. 
S tep 11 ; is oC ' and Mil? + 


let us say candidate 4 is 


Mil = 15. 

= 150. Therefore go to 


Step 12. 

Step. 12; INCMT = Mil? = 15, j Qut = v = 3 and 
m' ' = + INCMT = 150 

S o 

S tep 1 5s k is not equal to K at this step, therefore 
go to Step 14. 

Ste p.. 14; Set k = 1+1 = 2 and go to Step- 6. 



Cycle } from Step 6 to Step 14 , is repeated for next 
j , candidate which is V(2) = 8. 

° OU'G 

This yields I1TCMT =10 1. „ = 6 and Fl' =139. 

■ lor the last candidate V( 3 ) =9 the steps from 

6 to 11 are repeated. At Step 11 (MIH + ) happens to 

? 

he greater than M , therefore it goes to Step 13. 

Step 13 s k = 3 and K - 3, go to Step 15. 

Step 15: Mg - Mg = 153 - 139 = 14 

nicMf = 10 

IITCMT is less than (Kq, - li'), therefore go to 

to o 

Step 16. 

f 1 

Step 16 ; Make cyclic exchange jJ(8), 6 j on sepuence 
S = 3“5"l“Q"9-4“7“6-2. The new sequence will he 
S* = 3„5 -i»8-9-6-4-7-2. 

Step 17 2 Set S = S*= 3-5-l-8-3~S~4-7-2 and go to Step 2. 

This completes one iteration of the algorithm 
yielding a reduction of 4 units in the makespan. 

In the second iteration a further reduction of 3 units 
is realized. 

After the third iteration it stops without any 


more reduction. 
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So, finally we get the following sequence as 
the output of proposed heuristic. 

S = 3-5 -8-1-9 “6 -4 -7" 2 

Makespan = 146 

Total reduction = 153 -- 146 = 7 


CPU time 


0.15 seconds. 



CHAPTER Y 


PERFORMANCE 


VALUATION AND CONCLUSIONS 


5.1 Pe rforman ce Jjm^luatioii ; 

In this section s the results of experimental 
investigations carried out to evaluate the performance 
of the proposed heuristic procedure are presented. The 
evaluation is primarily done in terms of the number- of 
times the heuristic succeeds in improving the initial 
sequence and the amount of reduction in the makespan 
time. 

The experimental investigation was divided into 

two parts. firstly, an attempt was made to solve a 

standard problem available in the literature with known 

optimal solution. The 5-machine , 10- job flowshop scheduling 

* 

problem given in the book by Conway et.al. , was solved 
using the proposed heuristic. Conway, et al. , have 
claimed this to be the hardest flowshop scheduling 
problem available in the literature for solving through 
the use of heuristic procedures. The objective was to 

see as to how far the proposed heuristic brings the 

% 

initial suboptimal sequence close to the optimal solution. 
The initial suboptimal sequence was generated using 
P etro v* s algo r x thm._ __ 

* Conwajr, R, ¥. , Maxwell, W.L., and Miller, L.W., Theory of 
Scheduling-, Addison-Ueslcy, Reading, Mass. , 1967, p. 94. 
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The second part of the investigation involved the 
testing of the proposed heuristic on a set of randomly- 
generated problems of varied sizes. In all 13 problem 
sizes varying from 3-machine , 4™ .job to 7-machine, 10- job 
were considered. Each problem size was repeated five 
times. The processing times on the machines for various 
jobs were obtained using random number generators. Again, 
Petrov’s algorithm was used for the generation of initial 
s equenc e , 


5.1.1 E valuation of the Stand ard Prob lem s 

The processing ^ time matrix of the standard problem 
stated earlier is given below; 


T 


15 7 8 
2 9 6 9 
9 7 3 9 


5 7 


Ci 


6 



2 3.0 7 9 1 1 | 

3 4 7 4 3 1 j 


The optimum make span time for this problem has been reported 


as 6(5. 


The Petrov's heuristic gives the following 


sequence. 

S = 1-5-3-2-6-7-4-8-9-10 with a makespan of 69. 

Using Petrov’s solution as the initial sequence, 
the proposed heuristic in the first iteration gave a 
makespan value of 66. The corresponding sequence was 
l_2-5 -3-6 -7 -4-8-9-10. 



Thus for the 'hardest problem' the proposed 
heuristic succeeded in arriving at the optimal sequence 
ill just one iteration. The author has not been in a 
position to extensively test the abilitj/- of the nroposed 
heuristic to generate optimal sequence for want of large- 
sized standard problems with known optimal solutions. 

5.1.2 Eval Lu ation on Ra ndomly Gen er a te d Pro ble ms • 

The results obtained on average make span reduction 
for the randomly generated problems of various sizes are 
presented in Table 5.1= It was observed that the proposed 
heuristic improved the initial sequence for almost 70 per 
cent of the problems and the average improvement was 4 per 
cent. Purther, it was noted that the tendency to improve 
upon the initial solution increased with the increase in 
the problem size. 

In Table 5.1? the average computational requirements 
of the proposed heuristic for various problem sizes is 
also presented. The problems were solved on DEC 1090 
computer system and for each problem size a set of 5 problems 
were considered. It is observed that the computational 
requirements for the proposed heuristic are reasonable even 
for problems involving 7 machines and 10 jobs. 



5 • 2 Con e : lu s ions ; 

In this thesis, we have presented a heuristic 
methodology for improving upon an initial sequence for a 
static flow shop problem. The objective has been to exploit 
the inherent characteristics of the initial subopt imal 
sequence. The concepts of idle time, slack time, critical 
path and cyclic pairwise exchange have been utilized 
advantageously for the development of the proposed heuris- 
tic methodology. The experimental investigations have 
indicated that the heuristic algorithm improves the initial 
solution for most of the problems. The tendency to yield 
improved sequence increases with an increase in problem 
size, further, it is observed that the methodology is 
computationally efficient . 


Table 5.1 


SI. 

lo. 

Problem 

Size 

Average : 
make span 

reduction m 
time, p ere ent 

Average C 
( Secor 

1 . 

3x4 


0.03 

0.05 

2. 

3x5 


2.45 

0.08 

3. 

3x8 


2.90 

0.09 

4. 

4x5 


l P)£ 

0.07 

5. 

4x6 


3.30 

0.08 

6. 

4x7 


2. 83 

0.09 

7. 

5x7 


4.96 

0.11 

8. 

5x8 


4.02 

0.12 

9. 

5x3 


5.86 

0.15 

10. 

6x4 


4.17 

0.08 

11. 

6x3 


5.31 

0.20 

9 

C\i 

H 

7x9 


7.12 

0.21 

13. 

7x10 


6.13 

0.25 


?1 


tune 
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APPENDIX A 


PETROV'S ALGORITHM 


Given a processing time matrix T = 
m machine, n job flowshop problem. 

Define the corresponding 2-machine, 


fh.j] £o - 

n~job problems 



“1 

£ 

i=l 



the processing time for the 
j-th job on machine 1. 



m 

Y1 t. • = the processing time for the j-th 

i=m 2 ’ J 

job on machine 2. 


If ra is an even number then, 

= §r and m 0 = m. + 1 

If m is an odd number then 

m-j = and m2 = 


The following steps are applied to the above 
mentioned 2-machine problems 

/ 1 

St ep 1 ; Determine a column matrix | C^- > where 

°i = 1 - A 

Step.. 2 ; Apply rule No. 1; 

Arrange all the jobs, having Ch ^ 0, in increasing 

value of C-. s and the jobs, having C,< 0, in decreasing value 
j J 

of Ch. Call the sequence, hence generated, S-^ 

J 



Apply Rule Ho. II: 


Step, ji 

Arrange all the jobs In decreasing value of C- 

3 

generating a sequence S 2 . 

Step 4 • Select , from and 82 * the one which gives 
minimum make span. 

In addition to this, if there are some zero 
elements in the T‘ matrix then two more sequences are 
generated as follows: 


Define , 


T 

b 


”1 

rl. t. 
i-1 1 ’ J 


? 

3 . 


m 

X! ■ 

? . 
i=m. 


3 


- 0 

where j IT. is the number of machines out of first m^ machines 
having positive processing tines for job j and, IT - Is the 

J 

number of machines out of last ( n - m 2 +l) machines having 
positive processing times for job j. 

Apply .Rule Ho. 1 and 2 mentioned above to Aj and 

I 

r_. to generate two more sequences S* and S. . 

JO ■ 0 s- 

Seleot the best out of S-^, S 2 , S,. and S^. 



